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In the high density, low temperature limit, Quantum Chromodynamics exhibits a 
transition to phases characterized by color superconductivity and energy gaps in 
the fermion spectra. We review some fundamental results obtained in this area 
and in particular we describe the low energy effective lagrangian describing the 
motion of the quasi-particles in the high density medium (High Density Effective 
Theory). 



1 Lecture I: Color superconductivity 

1.1 Introduction 

At high baryonic densities and small temperatures the color interaction favors 
the formation of quark-quark condensates in the color attractive antisymmet- 
ric channel: 

A = (^ a C^)e a ^e^ k n ik jt (1) 

(a, j3, 7 = 1, 2, 3 color indices; i, j, k = 1, 2, 3 flavor indices). The condensates 
(|l|) depend on the matrix f2 and act as order parameters of new phases where 
the SU (3) c color symmetry is spontaneously broken. The densities are so high 
that these phenomena might probably occur only in the core of neutron stars. 
Since this mechanism is similar to electron superconductivity it is referred 
to as color superconductivity (CSC) □. This is one of the most fascinating 
advances in^Quantum Chromo Dynamics (QCD) in recent years (for 
reviews see □) . The aim of these lectures is to describe an approach ta JLhis 
aspect of QCD that is based on the method of the effective lagrangians &LrEltJ 
and to stress possible astrophysical consequences. 

Different phenomena take place depending on the value of the order pa- 
rameter ([j]). One could have: f2 7 fc = <5 7 3<5fc3, which corresponds to 

e^ey < VfaCV^ >= A<F 3 (2) 

where the sum over the flavor indices run from 1 to 2 and ip represents a left 
handed 2-component Weyl spinor (the right handed field satisfies a similar 
relation) ; moreover a sum over spinor indices is understood and C = ia% . This 
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case correspond to decoupling of the strange quark Qni* — * oo; m u — rrid = 0) 
and is called 2SC model. From dynamical analyses □ □ one knows that, for fx 
sufficiently large, the condensate ^ is non vanishing. Therefore it breaks the 
original symmetry group SU(3) c (g>SU(2) L <g l SU(2) R <giU(l)B down to S£/(2) c <g) 
SU (2) l ® SU (2) fl® Z2. The chiral group remains unbroken, while the original 
color symmetry group is broken to SU(2) C , with generators T A corresponding 
to the generators T 1 , T 2 ,T 3 of SU (3) c ; an unbroken C7 (1) ^ also remains. As a 
consequence, three gluons remain massless whereas the remaining five acquire 
a mass. One can construct an effective theory to describe the emergence of 
the unbroken subgroup SU (2) c and the low energy excitations, much in the 
same way as one builds up chiral effective lagrangian with effective fields at 
zero density. For the two flavor case this development can be found in HI . 

For the three flavor case [m u — m<j = m s — 0) the following case has 
been widely discussed a: 

3 

MaV&> = -<$M> = A E ^PKCijK ■ (3) 

K=l 

The condensate ||) breaks the original symmetry group SU (3) c CS> 5/7 (3) l ® 
SU(3) r <8)U(1)b downto SU(3) c+L+R <g>Z 2 . Both the chiral group, U{l) B and 
the color symmetry arc broken but a diagonal SU(3) subgroup remains unbro- 
ken in a way that locks together color and flavor (Color-Flavor-Locking=CFL 
model). There are 17 broken generators; since there is a broken gauge group, 8 
of these generators correspond to 8 longitudinal degrees of freedom of the glu- 
ons, because the gauge bosons acquire a mass; there are 9 Nambu-Godstone 
Bosons (NGB's) organized in an octet associated to the breaking of the flavor 
group and in a singlet associated to the breaking of the baryonic number—The 
effective theory describing the NGB for the CFL model was studied in til. 

Another interesting possibility arises when there is a difference Sfj, between 
the chemical potentials of the tjssip gapped quarks i-tfor references see 113; see 
also EiE-a and for earlier works As shown in 113 when the two fermions 
have different chemical potentials [i\ ^ /i2, for S/i of the order of the gap 
the vacuum state is characterized by a non vanishing expectation value of 
a quark bilinear which breaks translational and rotational invariance. The 
appearance of this condensate is a consequence of the fact that for fj,\ ^ [ii, 
and in a given range of <5/i = \fi% — (j,^, the formation of a Cooper pair with a 
total momentum pi + j>2 = 2<f ^ is energetically favored in comparison with 
the—normaLBCS state_^The possible form of these condensates is discussed 
in E3 and (see also E3); it is worthwhile to note that, for simplicity, these 
authors assume only two flavors. In E3 the ansatz of a plane wave behavior 
for the condensate is made : cx e 2lq x ; in 113 it is shown that the configuration 
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that is energetically favored is a face centered cubic crystal. For simplicity I 
will consider here only the plane wave: 

- < 0| eii e a ^ ia (£)C^(£)|0 >= 2Y h A e 2i * 3 ; (4) 

besides the scalar condensate (^) there is also a vector condensate that can be 
however neglected since it is numerically small. The 2SC and/or the LOFF 
phases might exist for intermediate values of the chemical potentials, while 
for very high /i the CFL phase should set in. 

The phase transition and the non vanishing condensates result from a 
mechanism analogous to the formation of an electron Cooper pair in a BCS 
superconductor. At T — the only QCD interactions are those involving 
fermions near the Fermi surface. Quarks inside the Fermi sphere cannot in- 
teract because of the Pauli principle, unless the interactions involve large 
momentum exchanges. In this way the quarks can escape the Fermi surface, 
but these processes are disfavored, as large momentum transfers imply small 
couplings due to the asymptotic freedom property of QCD. Even though in- 
teractions of fermions near the Fermi surface involve momenta of the order 
of fx, their effects are not necessarily negligible. As a matter of fact, even 
a small attractive interaction between fermions near the Fermi surface and 
carrying opposite momenta can create.an instability and give rise to coherent 
effects. This is what really happens tffi and the result is the formation of a 
diquark condensate, as expressed by @, (||) or (Q). We stress again that the 
only relevant fcrmion degrees of freedom are therefore those near the Fermi 
surface. In □ an effective theory for the CFL model was discussed, based on 
the approximation of the neglect of the negative energy states. This results 
in a rather terse formalism displaying as a characteristic note the existence of 
a Fermi velocity superselection rule and effective velocity-dependent fermion 
fields. We will refer to this effective lagrangian as the High Density Effective 
Theory (HDET). In 13 the 2SC model has been studied by the same formal- 
ism, while in □ this effective theory has been applied to the crystalline color 
superconducting phase E3, the so-called LOFF £3 phase. 

The aim of these lectures is to review some developments in the description 
of Color Super Conductivity that are based on the HDET approach. This will 
be mainly done in the second lecture. Here I wish to discuss some possible 
astrophysical implications of CSC, in particular in the LOFF phase. 

1.2 Astrophysical implications of the LOFF phase 

Besides its theoretical interest for the study of the phase structure of QCD 
theory, the crystalline phase may result relevant for astrophysical dense sys- 
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terns, in particular in the explanation of the glitches in the pulsars. 

The pulsars are rapidly rotating stellar objects, characterized by the pres- 
ence of strong magnetic fields and by an almost continuous conversion of ro- 
tational energy into electromagnetic radiation. The rotation periods can vary 
in the range 10~ 3 sec up to a few seconds; these periods increase slowly and 
never decrease except for occasional glitches, when the pulsar spins up with 
a variation in frequency that can be (5f2/f2 sa 10~ 6 or smaller. Glitches are a 
typical phenomenon of the pulsars, in the sense that probably all the pulsar 
have glitches. 

Pulsar are commonly identified with neutron stars; these compact stars 
are characterized by a rather complex structure comprising a core, an inter- 
mediate region with superfiuid neutrons and a metallic crust. The ordinary 
explanation of the glitches is based on the idea that these sudden jumps of 
the rotational frequency are due to the movements outwards of rotational vor- 
tices in the neutron superfiuid and their interaction with the crust. This is 
one of the main reasons that allow the identification of pulsars with neutron 
stars, as only neutron stars are supposed to have a metallic crust. Since the 
conventional models for glitches may be not familiar to an audience of nuclear 
and/or high energy physicists, I will briefly review them in the sequel. 

It is known that a boson liquid at T ~ forms a condensate whose 
wavefunction 

So = yfcfaPjeWrt (5) 

has a macroscopic meaning, due to the large number of particles in it. Also 
the probability current density 

J corid = ^(S VHS-3SVH ) = A noV $ (6) 
Ira m 

has such macroscopic meaning. Since j con d = nov s , where v s is the condensate 
velocity that can be identified with the superfiuid velocity, we get 

h - 

■ V$ . (7) 



rn 



The consequence of (m) is 



7 



• d£ = , (8) 



if the domain where the curve 7 lies is simply connected. Given the arbitrari- 
ness of 7, a different way to write this result is 

VAw>0. (9) 
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Let us now suppose that the vessel containing the liquid is put in rotation 
with angular velocity Q. A consequence of @ is that the superfluid cannot 
rotate; to the same result one arrives by noting that the vessel cannot commu- 
nicate the rotation to the superfluid component, as there is no friction between 
the recipient and the liquid. However it can be shown that the absence of ro- 
tation in the superfluid does not correspond to a state of minimal energy. As 
a matter of fact, if E and L are energy and angular momentum in an inertial 
frame, the energy as computed in the rotating frame is E rot — E — L ■ Q. If 
is sufficiently high, then a lower energy can be achieved with 

L-Q>0 (10) 

instead of L ■ f2 = that corresponds to the absence of rotation 111. 

Let us suppose for simplicity that the curve 7 lies in a plane. For the 
result (|l^) to be compatible with (||) one has to suppose that inside the curve 
7 there is one point where (|^) is violated. Physically this would correspond to 
the presence of a point with a normal, not superfluid component. Since now 
the domain is not simply connected, is substituted by 

v s ■ d£ = 2irnK , (11) 

'7 

where the integer n is a winding number counting the number of times the 
curve goes around the singular point, and n is a constant with dimensions of 
vorticity, i.e. \L} 2 ■ [T] _1 . To the same result one could arrive by noting that, 
by virtue of (0) , as <& and <f> + 2mr correspond to the same wavefunction, one 
may have 

I V$-df = 27rn, (12) 

J -V 

which shows that k — h/m; k is called quantum of vorticity. 

Clearly we can repeat the argument for any plane parallel to the previous 
one; we therefore conclude that there is an entire line (vortex line) of singular 
points. If this line is a straight line, then v s will be perpendicular to the vortex 
line and also perpendicular to the radius joining the singular point and the 
point at which we compute v s . At a distance r from the singular point one 
has 

v s = — , (13) 
r 

as can be immediately seen from (|ll[). More generally: 

d£AR 



(14) 
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where R is the distance vector from the vortex line (v.l.) to the point at which 
we compute the superfluid velocity. 

We note some properties of the vortex lines. First, the integral ( p"l| ) is 
independent of 7, provided the second curve contains the singular point; this 
result follows from the Stokes theorem. Second, the vortex line must be closed 
or it must stop at the boundary; were it open, for example at a point P, one 
might construct a surface £ lying on the contour 7 but large enough as to 
have no intersection with the vortex line; therefore by the Stokes theorem one 
would get that the vorticity constructed by v s would be zero. If the vortex 
line stops at the boundary and the vessel is rigid then the v.l. is pinned at 
the boundary, as it will be discussed in more detail below. 

Let us compute the critical angular velocity for the formation of the 
first vortex line. The formation of a vortex line changes the energy E rot by 

the amount AE rot = AE — A ( £ • Q ) . We have AE = / ^-dV. Here the 



2 

integral is over the entire volume of the vessel that, for simplicity, we assume 
to be a cylinder of height b and radius R; on the other hand v s is given by 
the expression valid for a vortex line, i.e. (|l3|). Therefore 

AE = —2tt v 2 rdr = bp s TT — — In — (15) 
2 J fci a 

where a is a cutoff of the order of the interatomic distances, at which the 
macroscopic description fails down. The minimal energy E rot is obtained for 

L parallel to f2, with L = / p s v s rdV = birR 2 — p s . Vortex lines appear if 

J m 

AE rot < 0, i.e. if 

— ln--r!i? 2 <0, (16) 
m a 

Ti R 

which corresponds to Q > Q, C rit = — t^t In — ■ Incidentally we note that, 

mR z a 

ceteris paribus, vortex lines with n = 1 are more stable than those with 



n > 1, as the positive term in (16) has an extra power of n. Therefore from 
now on we put n = 1. 

What happens when 3> f2 cr jt? Clearly we expect there will be several 
v.l.'s. During the rotation these vortex lines follow the rotational motion 
of the vessel, which is clear because they are pinned at the boundary of the 
superfluid^]. Their motion imitate the motion of the liquid as a whole, as it can 



a For rotations around an axis, the vortex lines are, by symmetry, straight lines parallel to 
the rotation axis. For motion inside holes, slits , etc. there can be closed v.l.'s that are 
called vortex rings. 
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be seen by the following argument: in the formula for AE cr n one can forget, 
for large fi, the first term AE in comparison to the second one. Therefore 
minimizing AE cr i t correspond to maximize the angular momentum L, which 
is obtained if the liquid moves as a whole. A consequence of this feature of 
the motion of the v.l.'s is that also for the superfluid one can use the well 

known hydrodynamical formula Q = —\7Av s , which strictly speaking is valid 

only for the normal component v n ; it can be used here only because it refers 
to the rotation of the superfluid as a whole. 

The number of vortex lines that are present in the superfluid is propor- 
tional to Q, according to the formula (iV=number of lines per unit area): 

which shows that with increasing N the v.l.'s tend to fill in all the space. 
To prove (0) we consider a large closed curve C encircling the area A and 
containing in its interior NA v.l.'s. One has 



= NA2itk (18) 



but also 



d£-v s = / da ■ V A v s = 2 AO, (19) 

C J A 

from which (0) follows. As an example we can evaluate TV for the pulsar 
in the Crab nebula. Here m = 2mjv (the condensate is formed by neutral 
bosons: pairs of neutrons) and £1 = Vl pu i sar gives N ~ 1.9 x 10 5 cm~ 2 with an 
average distance between vortex lines d ~ TV -1 / 2 ~ 10~ 2 cto. 

Let us consider again eqns. (0) and JlS}). If v(r) is the number of vortices 
per unit area at a distance r from the rotation axis, they give, if v = v s is the 
superfluid velocity, 

<fdi-v=[ dS ■ V A v = 2ttk f 2irr'v{r')dr' . (20) 
J Jo Jo 

We put 

k = 2ttk = — — (21) 



2m n 



and write (EOh as follows: 



2irr 2 n(r) = k [ 2nr'u(r')dr' , (22) 
Jo 
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which implies 



BCl 

kv(r) = 20(r) + r — . (23) 



Since the total number of v.l.'s is conserved, one has 

dv 



.., , V-K r )=0 (24) 
at 



where v r is the radial component of the superfluid velocity. We write ( P2| ) as 

2i:r 2 Cl{r) = k [ v dS (25) 



o 



and take the time derivative, using (|23|) to get 



2Trr 2 — = -k I dSV(vv r ) . (26) 







2 dn 

Using the Gauss theorem one gets 2nr — — = —k2nrvv r , i.e. 

at 

^ f 2m+r d -f]^. (27) 



dt r \ d 



r I r 



Eq. ( p7| ) shows that the only possibility for the superfluid to change its angular 
velocity (CI ^ 0) is by means of a radial motion, i.e. v r ^ 0. 

Let us now consider a rotating superfluid-in contact with rotating normal 
matter on which an external torque is acting E3. We denote by I c , Q c moment 
of inertia and angular velocity of the normal components that, in a neutron 
star, includes the crust and possibly other normal components. The equation 
of motion of the normal component is 

J c fi c (t) = M ext + M mt . (28) 

Besides the external torque M ext , basically related to the spin down of the 
pulsar (or the steady accretion in binary pulsars), we have included internal 
torque M mt : 



M int = ~ / dl p tt(r,t) (29) 



due to the interaction with the superfluid. Eqs.(|27j|29|) are the equations of 
motion for the angular velocities f2 and O c (superfluid and crust). The two 
velocities are coupled not only through Mint, but also by v r , because we will 
show below that v r = f(Cl — Cl c ). We note again that fundamental for this 
model is the existence of radial motion, for, if v r = 0, then CI = const, and 
only Cl c changes, due to the external torque alone. 
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In the neutron star, superfluid neutrons (in Cooper pairs) coexist with 
nuclei of the crust. Also in the crust there are superfluid neutrons, but they 
are characterized by a different (and smaller) A. For superfluid neutrons in 
the volume V, the total energy can be estimated as follows: 

V fPF + ^/E F dQ y A 2 k 3 

E ^wwL pdp - = ^^f (30) 

where pp = frkp; we integrate over half the solid angle as the superfluid 
neutrons only appear in pairs; we have taken into account that only neutrons 
in a shell of thickness A 2 /Ep participate in the pairing. 

E is also approximately given by the difference between the energies of 
superfluid neutrons outside the vortex line and neutrons inside, because those 
inside the vortex core have A — > 0. Therefore neutrons inside the volume 
V of the vortex core are repelled from going outside the vortex towards the 
superfluid phase as it would cost more energy. However, if neutron rich nuclei 
are present, the repulsion will be less important, as A c , the gap for superfluid 
neutrons in the nuclei, is much smaller than A s , the gap of superfluid neutrons; 
therefore there will be a force pulling the vortex toward the nuclei; the pinning 
energy per nucleus will be 



SE p 



V 



A 2 fc|A fA 2 k F 



TT Ep I \TT 2 Ep 



V ( A 2 k F 
7 8 V^E-p ] (31) 



with 7 ~ 1. 

Let now £ be the coherence length, i.e. the spatial extension of the Cooper 
pair, £ = — it can be proved that it gives an estimate of the radius of the 

7rA 

vortex core. The maximum pinning force will be obtained, if 2£ < b (b the 
average distance between the nuclei) when the vortex passes through one 
layer of the lattice; the average distance between vortex core neutrons and 
superfluid neutrons is of the order of £ and therefore the maximum force is 
SE 

F p ~ — and the maximum force per unit lenght of vortex line (b is also the 
average distance between two consecutive pinning centers) is 

U^- (32) 

Let us finally discuss a possible mechanism for the formation of glitches 0; 
we consider the rotating neutron star with superfluid neutrons in its interior 
and a metallic crust, which is a simplified model, but adequate to our purposes. 
We distinguish between two angular velocities: the superfluid velocity £1 and 
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the crust velocity fl c . Let us suppose that they are initially equal, which 
is a consequence of the pinning. Due to the spinning down of the star, fi c 
decreases; as long the vortex cores are pinned to the crust lattice, the neutron 
superfluid cannot spin down, because the radial motion is forbidden. There 
is therefore a relative velocity of the superfluid with respect to the pinned 
vortex core because Q > fl c : 

8v = (Q-Ci c )Ar. (33) 

The interaction between the normal matter in the core of the v.l. and the rest 
of normal matter (nuclei in the lattice, electrons, etc.) produces a Magnus 
force per unit length given by 

f = pkA5v, (34) 

where k is the quantum of vorticity and the direction of k coincides with the 
rotation axis. / is the force exerted on the vortex line; as it cannot be larger 
than f p there is a maximum difference of angular velocity that the system can 
maintain: 

^ = (^)^ = ^ = ^ (^) 

If u> < uj cr the vortices remain pinned at the lattice sites instead of flowing 
with the superfluid as they generally in superfluid (see discussion above) . On 
the contrary, if u) > oj cri the hydrodynamical forces arising from the mismatch 
between the two angular velocities ultimately break the crust and produce the 
conditions for the glitch. _ . 

We stop here this introduction to the standard model for glitches; see 113 
for more details. The relevance for CSC is that the LOFF phase provides 
a lattice structure independently of the crust. Even in quark stars, if one 
is in a LOFF phase, one has a crystal structure: a lattice characterized by 
a geometric array where the gap parameter varies periodically. This would 
avoid the frequently raised objection by which one excludes the existence of 
strange stars because, if the strange matter exists, quark stars should be rather 
common, in contrast with the widespread appearance of glitches in pulsars. 
Therefore, if the color crystalline structure is able to produce glitches, the 
argument in favor of the existence of strange stars would be reinforced. 

In a more conservative vein one can also imagine that the LOFF phase be 
realized in the inner core of a neutron star; in this case the crystalline color 
superconductivity could be partly responsible for the glitches of the pulsar. 
A detailed analysis of this scenario is however premature as one should first 
complete the study of the LOFF phase by including the third quark and by 
sorting out the exact form of the color lattice Il3'0. 
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2 Lecture II: High Density Effective Theory 

In this section we derive the general formalism of High Density Effective The- 
ory and we present an example of its use. This effective theorv-shows some 
resemblance with the Heavy Quark -Effective Theory (see e.g. t3); it is dis- 
cussed in detail in the review paper 0. 



2.1 General Formalism 

The main idea of the effective theory is the observation that the quarks par- 
ticipating in the dynamics have large (~ /i) momenta. Wherefore one can 
introduce velocity dependent fields by extracting the large part \xv of this 
momentum. One starts with the Fourier decomposition of the quark field 



= I 7^r^~ iv ' x m > (36) 



d A p 

and introduces the quark velocity by 

p" = + 1 , (37) 

where v» = (0, v) with \v\ = 1. Let us put l» = (£°, t) and £ = vt\\ + i± with 
t±_ = £—(£■ v)v. We can always choose the velocity parallel to p, so that 
4 = and 

d 4 p = [? ( dn f d£n f + d£ = 4tt^ 2 f ^- [ d£» f + d£ ■ (38) 



d 2 £e- u - x ^{£), 
(39) 



47T _ 

In this way the Fourier decomposition (|36|) takes the form 

where ipff(£) are velocity-dependent fields. One can write 

1>{x) = E e- ipv ' x [M*) + iM*)] , (40) 

V 

where ip± are velocity dependent fields corresponding to positive and energy 
solutions of the Dirac equation. 

Let us now define V* — (1, v), — (1, — v), 7^ = (7 , (v ■ 7) v), j± = 
7^ — 7y\ Using simple algebraic relations involving the gamma matrices □ one 
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obtains 

Cd = J2 U>liV-Di/> + +il; i _(2ii + iV-D)il>- +(^+i£) M 7^V-+h.c.)] ; (41) 

V 

is the covariant derivative: = + igA^. We note that here quark 
fields are evaluated at the same Fermi velocity; off-diagonal terms are sub- 
leading due to the Riemann-Lebesgue lemma, as they are cancelled by the 
rapid oscillations of the exponential factor in the \x — > oo limit. One may 
call this phenomenon Fermi velocity supers election rule, in analogy with the 
behaviour of QCD in the ttlq — ■> oo limit, where the corresponding effective 
theory, the Heavy Quark Effective Theory exhibits a similar phenomenon E3. 
By the same analogy we may refer to the present effective theory as High 
Density Effective Theory (HDET), 

We can get rid of the negative energy solutions by integrating out the ip- 
fields in the generating functional; in this way we get 

C D ^C q = Y, W+iV-dip + + ipUV-dip- , (42) 

V 

where now ij)± are both positive energy solutions with ip± = ip±z The con- 
struction described above is valid for any theory describing massless fermions 
at high density provided one excludes degrees of freedom very far from the 
Fermi surface. As discussed in the first lecture however, for small temperature 
and high density the fermions are likely to be gapped due to the phenomenon 
of the color superconductivity. We shall examine here the modification in the 
formalism for the LOFF model, with the condensate 

A(x) = A exp{2ig • x} . (43) 

The effect of the non vanishing vacuum expectation value can be taken into 
account by adding to the lagrangian the term: 

A rp 

C A = ~— exp{2iq- x} e a ^i^ ia {x)Ci)iji{x) - (L -► R) + h.c. . (44) 

In order to introduce velocity dependent positive energy fields iptf i; i a with 
flavor i we decompose both fermion momenta as in @ and we have: 

C A = ^2 exp{ix ■ a(vi, v j; q)}e tj e al33 'ip^ i ; z . la (x)Cip-^ :j . j p{x) 

Vi ,Vj 

— {L — » R)+h.c. , (45) 
where a(vi, Vj, q) = 2q — fj,iVi — (J-jVj. We define 

M= o ' °f 1= n ' ( 46 ) 
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and perform the y, —* oo limit on a smeared amplitude as follows 

lim exp{ix- a(v x , v 2 , q)} = — [ dr exp{if ■ d(v Xl v 2 , q)}, (47) 

^oc V J v{3) 

where V(x) is a small volume centered at x. We evaluate ( f47| ) by taking q 
along the z — axis, so that we get: Vi ~ — Vj = —v by the x and y integrations, 
while the z-integration gives (ft = q/q) 



n :e l2 « hz 5 R [h(v-n)]^^S R [h(v-n)}. 



H • R- - ^ 

We have put R — q\Al\ where A£ is a smearing distance along the direction 

7T 

of q (\A£\ ~ — ). We have introduced the "fat delta" Sr(x) defined by 



Sr(x) 



sin(i?a;) 



(49) 



which, for large R, gives S R (x) — > S(x). Moreover in the \i — > oo limit h(x) 
5[i 



Therefore one has 



qx 



A 7T 

~2 R 



5 a [h(v-n)] J~] eije a p^. ia (x)Cil}_ &,jp{x) - (L -> i?) +h.. 



In an appropriate basis the effective lagrangian is 



v A=0 



A | ^ ic>ab V ■ d A A b 

A A b i 5ab V ■ d 



X 



(50) 



(51) 



where the matrix Aab is as follows: Aab — (Aor B = 4 or 5), and, for 
A,B = 0, ...,3: 

/l \ 
0-10 
0-10 
\0 0-1/ 



MB 



(52) 



AB 



having put 



Ae// = ~^ dR ^" >i ^ ■ 



From these equations one can derive the propagator for gapped quarks: 

1 



D AB {£) 



V ■£ V ■£ 



A 2 



V-i8 AB -Aab 
-Aab V-ISab 



(53) 



(54) 
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2.2 NGB and their parameters 

Both in the CFL model and in the LOFF model with two flavors there are 
Nambu-Goldstone Bosons, associated in one case to the breaking of internal 
symmetries and in the other to the breaking of space-symmetries (phonon). 
In order to derive an effective low energy lagrangian for the NGB's one can 
use the gradient expansion, where the NGB's are introduced as external fields 
and acquire a kinetic term, thus becoming dynamical fields, by integrating 
out the fermion fields. I shall describe in some detail the calculation for the 
phonon field in the LOFF phase (other examples can be found in El) . 

Ca in ( fjif ) explicitly breaks rotations and translations and induces a 
lattice structure given by plane waves. The crystal can fluctuate and its local 
deformations define one phonon field (j> that is the Nambu-Goldstone boson 
associated to the breaking of the translational symmetry. It is introduced by 
the substitution in (|43|) 



z + 



(55) 



with ((f>)o = 0. We are interested in an effective description of the field <\> in 
the low energy limit, i.e. for wavelengths much longer than the lattice spacing 
~ 1/q. In this limit the field <fi varies almost continuously and we can get rid 
of the lattice structure and use in the sequel the continuous notation. 
At the first order one gets the couplings 

C<t>^ = - ^ S R [h(v ■ n)} 1 -^ ejje^Tp^a^C -ipj^-y - (L —> R) + h.c. 



C<t>W,$ = £ S R [h(v ■ n)} — eije^ipi^a^Cipj.p-v - (L -> R) 



In the basis of the x fields one has: 



(56) 

h.c. 

(57) 



Here 



7rA 



v A=0 



ot 



— 5 R [h{v ■ n)]— — 6 R [h(v ■ ri]—— 



R 



f 



R 



/l \ 
0-10 
0-10 

\o 0-1/ 



(58) 



(59) 



AB 
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To compute the effective lagrangian for the phonon field we use the prop- 
agator given in Eq. fl54| ) and the interaction vertices in (|5|). The result of 
the calculation at the second order in the momentum expansion is 

^ w - « is£ E 5 0£ w - ft " J1(1 *■>)'/ dH 2Al wm"' P ■ 

(60) 

One can handle the fat delta according to the Fermi trick in the Golden Rule; 
in expressions involving the gap parameters one makes in the numerator the 
substitution 8a[h(x)] — > S[h(x)}; one fat delta is substituted by the Dirac delta 

while the other fat delta gives -^j — > \ yy e finally get the 

R R 

effective lagrangian in the form 



9J9 > o < v ■ n 



W/W = - t^tj > J\ v-n - — ) VyVvV^vA ■ (61) 



Here kn is kinematical factor of the order of 1 induced by the approximation 
of the Riemann-Lebesgue lemma u. The integration over fermi velocities can 
be easily performed and one obtains the effective lagrangian in the form 

C^,d^) = \(^-v\\V^\ 2 -vl\^^\ 2 ) , (62) 



(63, 

Here Vii0 = n(H ■ V)<p , Vj_0 = W(j> — Vm0. Moreover 

«i = J sin 2 6 q , (64) 



'|| — Uq 

and 

2 



v\ = cos 2 (9 g (65) 



cos 2 9 q = ( 5 -^\ . (66) 



V i 

Therefore the dispersion law for the phonon is anisotropic 



b One can easily control that the Goldstone theorem is satisfied and the phonon is massless. 
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